Joint Design of a Simultaneous Reflection and Transmission RIS in Mode-Switching Mode to Assist NOMA Systems

Simultaneous transmitting and reflecting reconfigurable intelligent surfaces (STAR-RISs) can reflect signals and transmissive signals simultaneously and can extend the coverage of signals. A conventional RIS mainly focuses on the case where the signal source and the target are on the same side. In this paper, a STAR-RIS-assisted non-orthogonal multiple access (NOMA) downlink communication system is considered to maximize the achievable rate for users by jointly optimizing the power-allocation coefficients, active beamforming, and STAR-RIS beamforming under the mode-switching (MS) protocol. The critical information of the channel is first extracted using the Uniform Manifold Approximation and Projection (UMAP) method. Based on the key extracted channel features, STAR-RIS elements and users are clustered individually using the fuzzy C-mean clustering (FCM) method. The alternating optimization method decomposes the original optimization problem into three sub-optimization problems. Finally, the sub-problems are converted to unconstrained optimization methods using penalty functions for the solution. Simulation results show that when the number of elements of RIS is 60, the achievable rate of the STAR-RIS–NOMA system is about 18% higher than that of the RIS–NOMA system.


Introduction
With the rapid development of modern mobile communication technologies, the capacity of communication networks will grow significantly in the next decade, and the high complexity of networks, the high cost of hardware, and the increasing consumption of resources have become the key issues facing wireless communication in the future [1]. Reconfigurable intelligent surfaces (RISs) are low-cost, programmable, easy to deploy, and low-energy, which makes RISs stand out in competition with other technologies. RISs and their variants [2,3] have become promising technologies for sixth-generation (6G) wireless networks. Simultaneous transmitting and reflecting (STAR)-RIS is a new technique in RISs; whereas a conventional RIS focuses on the case where the signal source and the target are located on the same side of the RIS [4][5][6], a STAR-RIS is capable of simultaneous transmission and reflection, and the target can be distributed on both sides of the signal source, allowing for full spatial coverage and enabling better-aided system communication [7][8][9].
Due to the significant capacity growth of future communication networks, nonorthogonal multiple access (NOMA) [10][11][12][13][14] has also received considerable attention as a promising technology to enhance spectrum efficiency and support large-scale connectivity. NOMA transmits the information of multiple users on the same sub-channel. Interference information is introduced, in contrast to orthogonal multiple access (OMA) systems, where each user has a separate resource in the time, frequency, or code domain [15]. NOMA works through sharing time-frequency resources by separating users into other domains [16].

System Mode
As shown in Figure 1, consider the downlink transmission of the STAR-RIS-NOMA system, where the base station (BS) communicates with M single-antenna users, the set of users is M {1, · · · , M}, and the direct link between the users and the BS is blocked. Assume that the BS is equipped with N t transmit antennas, while the STAR-RIS is equipped with N ris elements and the set of elements is N ris {1, · · · , N ris }. All perfect CSIs are available at the BS. H ∈ C N ris ×N t is the channel from the BS to the STAR-RIS, and g m ∈ C N ris ×1 is the channel from the STAR-RIS to the m-th user. We use fuzzy C-mean (FCM) clustering to group M users into L clusters. The set of users in the l-th cluster is denoted by L l . Thus, the number of users in the l-th cluster can be denoted by |L l |, where ∑ L l=1 |L l | = M. The detailed process on user clustering will be discussed in Section 3. This paper considers the user-achievable sum-rate maximization problems under the mode-switching (MS) protocol. Under the MS protocol, STAR-RIS elements are divided into two groups: one working in reflection mode (R mode) and the other in transmission mode (T mode). Let u t = where β t n ris , β r n ris ∈ {0, 1},β t n ris + β r n ris = 1 and φ t n ris , φ r n ris ∈ [0, 2π), ∀n ris ∈ N ris . The corresponding STAR-RIS transmission and reflection coefficient matrices are Φ t = diag u t and Φ r = diag(u r )., respectively.
In the STAR-RIS-NOMA system, v l denotes the active beamforming vector of cluster l, and s l,c denotes the signal sent by the BS to the c-th user of the l-th cluster. Therefore, the received signal of the c-th user of the l-th user cluster is given by y l,c = g H l,c Φ l,c Hv l √ η l,c s l,c + ∑ i =c,i∈L l g H l,c Φ l,c Hv l √ η l,i s l,i + ∑ j =l,j∈L ∑ i∈L j g H l,c Φ l,c Hv j √ η j,i s j,i + n l,c , where η l,c denotes the power-allocation coefficient of the c-th user of the l-th user -cluster, n l,c represents the additive Gaussian white noise, and n l,c ∼ CN 0, σ 2 , σ 2 is the noise power. Due to the presence of multiple users, each user uses the SIC sequentially to eliminate interference. This paper considers the user-achievable sum-rate maximization problem in the case of fixed decoding order. For simplicity and without loss of generality, we assume that the equivalent channel gain order in the l-th user cluster is Thus, the achievable rate for the c-th user of the l-th user cluster can be given as In the STAR-RIS-NOMA system, l v denotes the active beamforming vector of cluster l , and , lc s denotes the signal sent by the BS to the -th c user of the -th l cluster.
Therefore, the received signal of the -th c user of the -th l user cluster is given by

Problem Formulation
To improve the overall data rate, the joint power-allocation coefficients, active beamforming, and transmission and reflection beamforming optimization problems at the RIS are proposed to maximize the achievable rate for users. The problem formulations are β o n ris ∈ {0, 1}, β t n ris + β r n ris = 1, where o ∈ {t, r}. Constraint (4b) guarantees a minimum quality of service (QoS) requirement for each user, and constraint (4c) indicates a total transmit power budget of P max . Since the objective function is nonconvex, the constraint (4b)-(4d) is also nonconvex. Therefore, problem (4) is a nonconvex optimization problem, which is more difficult to solve directly. Compared with reflection-only RIS, the STAR-RIS has a more complex optimization problem because of the presence of both reflection and transmission modes.

Channel Feature Extraction and Clustering of Users and STAR-RIS Elements
The number of transmitting antennas arranged at the BS is usually high. The number of STAR-RIS elements is significant, resulting in a high dimensionality of the relevant channels, which is challenging to optimize directly. The Uniform Manifold Approximation and Projection (UMAP) [39] can identify the critical structures in the high-dimensional channel space and preserve these structures in the low-dimensional embedding. Therefore, we first process the channel from the BS to the STAR-RIS and the STAR-RIS to the user using UMAP to extract the critical information and then perform clustering for STAR-RIS elements and system users, respectively.

Channel Feature Extraction Using UMAP
According to the three assumptions of UMAP about the data, a manifold with a fuzzy topology can be modeled. The embedding is found by searching the low-dimensional projection of the data with the closest equivalent fuzzy topology. The process of channel feature extraction using UMAP is shown in Figure 2. The STAR-RIS-user channel matrix is To extract the features of the STAR-RIS-user channel, the first step is to establish the relationship between the highdimensional samples. In the STAR-RIS-user channel matrix G , the probability-density function between each point can be given as  . The relationship between the low-dimensional points can be established by Equation (7).
Establishing the relationship between a point and its P nearest neighbors in a high-dimensional space Establishing the relationship between a point and its P nearest neighbors in a low-dimensional space Optimal low-dimensional results are obtained by optimizing the crossentropy function The STAR-RIS-user channel matrix is G = [g 1 , . . . , g M ]. To extract the features of the STAR-RIS-user channel, the first step is to establish the relationship between the highdimensional samples. In the STAR-RIS-user channel matrix G, the probability-density function between each point can be given as where d(·) denotes the Euclidean distance between sample points and ρ (G) i is the distance of the nearest point to point i. To ensure the symmetry of the distance, the joint probability is established as p Fix the R (G) nearest neighbors of the point i. Solve for θ The relationship between the low-dimensional points can be established by Equation (7).
where a (G) , b (G) is the hyperparameter and g i is the downscaled channel. In order to describe the relationship between the high-dimensional sample points and the reduceddimensional sample points, the cross-entropy-loss function is constructed as (8) where g ij denotes the distance between the high-dimensional sample points and g ij denotes the distance between the low-dimensional sample points. Therefore, the low-dimensional space g i , corresponding to g i in the high-dimensional space, can be obtained by solving the following problem: min CE(g ij , g ij ) We can use (10) to solve problem (9) using a stochastic gradient descent (SGD) update to obtain the reduced channel matrix G = [ g 1 , . . . , g M ].
where α 1 is the step size. The channel from the BS to the STAR-RIS is H = h 1 , . . . , h N ris H . First, it is necessary to establish the relationship between the high-dimensional samples. The probability-density function between each point in the BS to the STAR-RIS channel matrix H can be given as where d(·) denotes the Euclidean distance between sample points, and ρ (H) i is fixed as the distance of the nearest point to point i. To ensure the symmetry of the distance, the joint probability is established: Fix the R (H) nearest neighbors of the point i. Solve for θ The relationship between the low-dimensional points can be established by (13).
where a (H) , b (H) is the hyperparameter and h i is the downscaled channel. To describe the relationship between the high-dimensional sample points and the downscaled lowdimensional sample points, a cross-entropy-loss function is constructed: where h ij denotes the distance between the high-dimensional sample points and h ij denotes the distance between the low-dimensional sample points. Therefore, h i in the high- dimensional space corresponds to the low-dimensional space h i , which can be obtained by solving the following problem: Problem (15) is solved by SGD using (16) to obtain the channel matrix H = h 1 , . . . , h N ris T of the reduced BS to the STAR-RIS.
where α 2 is the step size.

User Clustering
After extracting the STAR-RIS-user channel features using UMAP and then using FCM-clustering to subdivide the user clusters, it is possible to classify them more accurately with less complexity. The process of clustering users by FCM is shown in Figure 3. Using FCM, we divide the users into L clusters and the sub-cluster set into L = {1, . . . , L}. The users are thus classified by the optimization problem (17).
where χ (user) is an affiliation factor showing the importance of whether a sample belongs to a class or not, Z is the affiliation matrix z indicating the affiliation of sample i to class j, and K = [k 1 , . . . , k L ] is the cluster center. The constraint (17b) is combined into the objective function using the Lagrange multiplier method: where ϑ (user) m is the Lagrange multiplier. Taking the derivative of (18) with respect to z ji and k j , respectively, yields (19) and (20).
Solving Equations (19) and (20), we can obtain Equations (21) and (22) are iterated and, finally, the users are divided into L clusters. The set of users in the l-th user cluster is L l , and the number of users is |L l |.
Clustering Center User 0.5 Affiliation 0.6 0.7 Step 1: Initialize clustering centers and affiliation 0.5 Step 2: Update the clustering center 0.4 Step 3: Update affiliation 0.6 Loop steps 2 and 3 until convergence

STAR-RIS Element Clustering
. , x N ris ; then, also using FCM, divide the elements of the STAR-RIS into |γ| clusters and the set of sub-clusters into γ = {1, . . . , |γ|}. The process of FCM-clustering of STAR-RIS elements is shown in Figure 4, and the classification results of STAR-RIS elements are obtained by the optimization problem (23): where χ (ris) is an affiliation factor showing the importance of whether a sample belongs to a class or not, A is the affiliation matrix, a indicates the affiliation of sample i belonging to class j, and D = d 1 , . . . , d L ris is the cluster center. The constraint (23b) is combined into the objective function using the Lagrange multiplier method: where ϑ (ris) n is the Lagrange multiplier. Equation (24) can be obtained by deriving Equations (25) and (26) for a ji and d j respectively. Solving Equations (25) and (26), we obtain Equations (27) and (28) are iterated and, finally, the N ris elements of the STAR-RIS are divided into |γ| clusters. The set of elements in the r-th RIS component cluster is γ r , and the number of elements is |γ r |. Step 1: Initialize clustering centers and affiliation Step 2: Update the clustering center Step 3: Update affiliation Loop steps 2 and 3 until convergence  After the channel feature extraction and the clustering operation, the received of the -th c user of the -th l user cluster transmitted/reflected through the -th r RIS element cluster is We assume that the power-allocation coefficients satisfy . Thu achievable rate of the -th c user of the -th l user cluster transmitted/reflected th the -th r STAR-RIS cluster elements can be given as Equation (30).  After the channel feature extraction and the clustering operation, the received signal of the c-th user of the l-th user cluster transmitted/reflected through the r-th STAR-RIS element cluster is y r,l,c = g H r,l,c Φ r,l,c H r v r,l √ η r,l,c s r,l,c + ∑ i =c,i∈L l g H r,l,c Φ r,l,c H r v r,l √ η r,l,i s r,l,i We assume that the power-allocation coefficients satisfy ∑ c∈L l η r,l,c = 1. Thus, the achievable rate of the c-th user of the l-th user cluster transmitted/reflected through the r-th STAR-RIS cluster elements can be given as Equation (30).
After clustering the STAR-RIS elements and users separately, the STAR-RIS element clusters correspond to the user clusters, in which case β t n ris and β r n ris are determined. Therefore, to maximize the achievable sum rate of M users, problem (4) is reformulated as where o ∈ {t, r}. Let U r.l.c = u r,l,c u H r,l,c ,b r,l,c = diag g H r,l,c H r , B r.l.c = b r,l,c b H r,l,c and V r,l = v r,l v H r,l , then, the achievable rate is Problem (31) can be reformulated as Rank(U r,l,c ) = 1, (33e) where f 1 (η r,l,c , U r,l,c , V r,l ) is as in (34), f 2 (η r,l,c , U r,l,c , V r,l ) is as in Equation (35), and f 3 (η r,l,c , U r,l,c , V r,l ) is as in Equation (36). Problem (33) is still nonconvex, and we use SDR to relax the rank-one constraint and then use a first-order Taylor expansion to approximate the nonconvex part of the objective function with the constraint present.
−η r,l,c Tr(U r,l,c B r.l.c V r,l )

Solution of the Problem
In this section, we propose a new alternating optimization iterative algorithm to solve the joint optimization problem of power-allocation coefficients, active beamforming, and STAR-RIS passive beamforming.

Power-Allocation Coefficient Optimization
When fixing the STAR-RIS coefficients and the BS beamforming vectors, problem (33) can be reformulated as At this point, the constraint is convex, and the first-order Taylor expansion of f 1 (η r,l,c ) is where η (τ) r,l,c denotes the value of η r,l,c at the τ-th iteration and, according to the SCA principle, Equation (37a) can be replaced by Using the penalty-function method, problem (37) can be transformed into an unconstrained optimization problem: where ζ 1 is the penalty factor. Thus, the gradient-descent method can be used to solve for η r,l,c by using Equation (41).
where α 3 is the step size.
Using the SDR relaxation rank-one constraint (33d), the first-order Taylor expansion of f 1 (V r,l ) is where V (τ) r,l denotes the value of V r,l at the τ-th iteration. Thus, the objective function (42a) can be converted to After removing the constraints (42b) and (33c) using the penalty function, problem (42) can be converted to Problem (45) is a semidefinite programming problem (SDP), which is solved by using the BCD [40] method. First, problem (45) is converted into an unconstrained optimization problem: where 1 > 0 is the barrier factor. The key to applying the coordinate-descent method to the problem of (46) is the following block determinant: where V r,l = S x x T y , and therefore, problem (46) can be transformed by solving where f 1 x is (49), f 2 x is (50), and f 3 x is Equation (51). (49) The gradient-descent method is used to solve for x according to Equation (52).
where α 4 is the step size.

STAR-RIS Coefficient Optimization
Fixing the power-allocation coefficients and BS beamforming vectors, problem (33) can be formulated as Regarding the relaxation of constraint (33e) using SDR, the first-order Taylor expansion of f 1 (U r,l,c ) is where U where ζ 3 is the penalty factor. Problem (56) is an SDP problem, and the solution using the BCD method can convert (56) into an unconstrained problem.
where 2 > 0 is the barrier factor. The key to applying the coordinate-descent method to the problem of (57) is the following block determinant: (60) (61) The gradient-descent method is used to solve forx according to Equation (63).
where α 5 is the step size factor. In Algorithm 1, the complexity of computing the power allocation is

Numerical Results
In this section, numerical simulations are performed to evaluate the performance of the proposed algorithm. The specific simulation scenario is shown in Figure 5. Without loss of generality, it is assumed that the BS and the STAR-RIS are located at (0, 0, 50) and (0, 40, 50), respectively. The STAR-RIS-NOMA system has M = 12 users. The users are divided into L = 4 clusters with cluster center coordinates of (0, 20, 0), (0, 30, 15), (0, 50, 0), and (0, 60, −15), respectively, and the users are randomly distributed in clusters corresponding to a radius of 6. It is assumed that STAR-RIS is equipped with N ris = N h RIS N v ris elements, where N h RIS = 5 denotes the number of elements along the horizontal direction, N v RIS = 5 represents the number of elements along the vertical path, and STAR-RIS elements are divided into |γ| = 4 clusters. The narrowband quasi-static fading channel from the BS to the STAR-RIS and from the STAR-RIS to the user is modeled as a Rician-fading channel as follows: where ξ BR and ξ RU,l,c denote the distance between the BS and the STAR-RIS and the distance between the STAR-RIS and the user, respectively; ς o denotes the path loss at a reference distance of 1 m; ω BR and ω RU denote the corresponding path loss index; B BR and B RU denote the Rician factors; and H Los and g Los l,c are the determined line-of-sight (LoS) components. H NLos and g NLos l,c are the random non-line-of-sight (NLoS) components modeled using Rician fading. In order not to lose generality, assume that users have the same QoS requirements and set R min l,c = 0.1 bits/s/Hz . The system parameters used are shown in Table 1.
where  BR and  ,, RU l c denote the distance between the BS and the STAR-RIS and the distance between the STAR-RIS and the user, respectively;  o denotes the path loss at a ref-  Table 1.     Initialize penalty factor 10 −4 τ max Maximum number of inner-layer iterations 25 Figure 6 investigates the convergence of the power-allocation algorithm. With P max = 35 dBm, N ris = 25, and N t = 5, the user-achievable sum rate of the algorithm converges in about 11 iterations as the number of iterations increases.  Figure 7 shows the relationship between the user-achievable sum rate and the number of RIS elements. A total of four cases are considered. The first case is the use of the STAR-RIS in the NOMA system, where the number of elements of STAR-RIS is = ris NN . In the second case, two conventional RISs are used in the system: one for reflection and Figure 6. Convergence of the power-allocation algorithm. Figure 7 shows the relationship between the user-achievable sum rate and the number of RIS elements. A total of four cases are considered. The first case is the use of the STAR-RIS in the NOMA system, where the number of elements of STAR-RIS is N ris = N. In the second case, two conventional RISs are used in the system: one for reflection and the other for transmission, with N elements per RIS. The third case also uses conventional RISs, but each RIS has N 2 elements. The fourth type is the NOMA system without RIS assistance. The results show that the performance of the STAR-RIS is higher than that of the conventional RIS for the same number of RIS elements. The higher version of the dual RIS than the STAR-RIS is because the dual RIS has N reflected and N transmitted elements at the same moment, while the transmitted/reflected elements of the STAR-RIS add up to a total of N. Figure 6. Convergence of the power-allocation algorithm. Figure 7 shows the relationship between the user-achievable sum rate and the number of RIS elements. A total of four cases are considered. The first case is the use of the STAR-RIS in the NOMA system, where the number of elements of STAR-RIS is = ris NN . In the second case, two conventional RISs are used in the system: one for reflection and the other for transmission, with N elements per RIS. The third case also uses conventional RISs, but each RIS has 2 N elements. The fourth type is the NOMA system without RIS assistance. The results show that the performance of the STAR-RIS is higher than that of the conventional RIS for the same number of RIS elements. The higher version of the dual RIS than the STAR-RIS is because the dual RIS has N reflected and N transmitted elements at the same moment, while the transmitted/reflected elements of the STAR-RIS add up to a total of N .   Figure 8 shows the user-achievable sum rate versus the BS transmit antenna in the case of conventional RIS-assisted NOMA and STAR-RIS-assisted NOMA communications. The experimental results show that the use of STAR-RIS-assisted communication is superior to a conventional RIS when the number of BS antennas is the same, given the same number of RIS elements. It can be observed that using RIS in NOMA can significantly improve the system's performance. Furthermore, as the number of antennas at the BS increases, the performance of the RIS-NOMA system also increases. superior to a conventional RIS when the number of BS antennas is the same, given the same number of RIS elements. It can be observed that using RIS in NOMA can significantly improve the system's performance. Furthermore, as the number of antennas at the BS increases, the performance of the RIS-NOMA system also increases.  Figure 9 shows the user-achievable sum rate versus the total transmitted power Pmax during conventional RIS and STAR-RIS-assisted NOMA communications. The userachievable sum rate increases with the total transmission power. RIS-assisted communi-  Figure 9 shows the user-achievable sum rate versus the total transmitted power Pmax during conventional RIS and STAR-RIS-assisted NOMA communications. The user-achievable sum rate increases with the total transmission power. RIS-assisted communications significantly improve user achievability and rates over NOMA systems without RIS assistance. It was observed that the performance difference between the system using two conventional RISs and the STAR-RIS system was relatively small.  Figure 9 shows the user-achievable sum rate versus the total transmitted power Pmax during conventional RIS and STAR-RIS-assisted NOMA communications. The userachievable sum rate increases with the total transmission power. RIS-assisted communications significantly improve user achievability and rates over NOMA systems without RIS assistance. It was observed that the performance difference between the system using two conventional RISs and the STAR-RIS system was relatively small.  Figure 10 considers the relationship between the user-achievable sum rate and the minimum user rate for three cases: a STAR-RIS system with N elements, a RIS system with 2N elements, and a RIS system with N elements. The results show that the user- Figure 9. Achievable sum rate relative to total transmitting power for N ris = 25 and N t = 5. Figure 10 considers the relationship between the user-achievable sum rate and the minimum user rate for three cases: a STAR-RIS system with N elements, a RIS system with 2N elements, and a RIS system with N elements. The results show that the user-achievable sum rate decreases as the minimum user rate increases. In all three cases, the RIS system with 2N components has the highest user-achievable sum rate, and the STAR-RIS system with N elements performs better than the RIS system with N elements. achievable sum rate decreases as the minimum user rate increases. In all three cases, the RIS system with 2N components has the highest user-achievable sum rate, and the STAR-RIS system with N elements performs better than the RIS system with N elements.  Figure 11 shows how the proposed method in this paper compares with four other methods or schemes. The first scheme is where no RIS is configured in the system. The second scheme is where each cell of the STAR-RIS is randomly configured with a phase. The third scheme uses a pairwise-rising method [41] to determine the phase shift of the STAR-RIS. The fourth scheme is to use the Riemann conjugate-gradient method [42] to solve the phase-shift matrix of the STAR-RIS. As can be seen from the results in Figure 11, when the number of RIS elements is 60, the user-achievable sum rate of the proposed Figure 10. Achievable sum rate versus minimum user rate for P max = 35 dBm, N ris = 25 and N t = 5. Figure 11 shows how the proposed method in this paper compares with four other methods or schemes. The first scheme is where no RIS is configured in the system. The second scheme is where each cell of the STAR-RIS is randomly configured with a phase. The third scheme uses a pairwise-rising method [41] to determine the phase shift of the STAR-RIS. The fourth scheme is to use the Riemann conjugate-gradient method [42] to solve the phase-shift matrix of the STAR-RIS. As can be seen from the results in Figure 11, when the number of RIS elements is 60, the user-achievable sum rate of the proposed method in this paper is 13.1 bits/s/Hz. The achievable sum rate of the random phase-shift method is only 10.8 bits/s/Hz. The achievable sum rates using the pairwise-ascent method and the Riemann conjugate-gradient method are 12.8 bits/s/Hz and 12 bits/s/Hz, respectively.  Figure 11 shows how the proposed method in this paper compares with four other methods or schemes. The first scheme is where no RIS is configured in the system. The second scheme is where each cell of the STAR-RIS is randomly configured with a phase. The third scheme uses a pairwise-rising method [41] to determine the phase shift of the STAR-RIS. The fourth scheme is to use the Riemann conjugate-gradient method [42] to solve the phase-shift matrix of the STAR-RIS. As can be seen from the results in Figure 11, when the number of RIS elements is 60, the user-achievable sum rate of the proposed method in this paper is 13.1 bits/s/Hz. The achievable sum rate of the random phase-shift method is only 10.8 bits/s/Hz. The achievable sum rates using the pairwise-ascent method and the Riemann conjugate-gradient method are 12.8 bits/s/Hz and 12 bits/s/Hz, respectively. Figure 11. Achievable sum rate relative to the number of RIS elements under different methods. Figure 11. Achievable sum rate relative to the number of RIS elements under different methods. Figure 12 shows the variation of achievable users and rates with the number of antennas at the BS. When the number of BS transmitting antennas is 10, the user-achievable sum rate of the random phase shift is 7.3 bits/s/Hz. The user-achievable sum rates of the pairwise-ascent method and Riemann conjugate-gradient method are 10.8 bits/s/Hz and 11.5 bits/s/Hz, respectively. The achievable sum rate of the proposed method is 12.6 bits/s/Hz. Among the five schemes, the worst performance is achieved without RIS assistance, followed by the random RIS phase shift. The Riemann conjugate-gradient method is closer to the case of the pairwise-ascent method, and the proposed method in this paper has the best performance. Therefore, this shows the effectiveness of the proposed algorithm.
Sensors 2023, 23, x FOR PEER REVIEW 21 of 23 Figure 12 shows the variation of achievable users and rates with the number of antennas at the BS. When the number of BS transmitting antennas is 10, the user-achievable sum rate of the random phase shift is 7.3 bits/s/Hz. The user-achievable sum rates of the pairwise-ascent method and Riemann conjugate-gradient method are 10.8 bits/s/Hz and 11.5 bits/s/Hz, respectively. The achievable sum rate of the proposed method is 12.6 bits/s/Hz. Among the five schemes, the worst performance is achieved without RIS assistance, followed by the random RIS phase shift. The Riemann conjugate-gradient method is closer to the case of the pairwise-ascent method, and the proposed method in this paper has the best performance. Therefore, this shows the effectiveness of the proposed algorithm.

Conclusions
In this paper, we study the downlink transmission system of the STAR-RIS-NOMA. A new alternate optimization method is proposed for the MS protocol to jointly optimize the power-allocation coefficient of the user, active beamforming at the BS, and passive

Conclusions
In this paper, we study the downlink transmission system of the STAR-RIS-NOMA. A new alternate optimization method is proposed for the MS protocol to jointly optimize the power-allocation coefficient of the user, active beamforming at the BS, and passive beamforming of the STAR-RIS so that the achievable rate of the user is maximized. The channel is first feature-extracted, and the user and STAR-RIS elements are clustered according to the extracted features to reduce the solution difficulty. The SDR technique and penalty function are used to solve the nonconvex factors present in the problem. The effectiveness of the proposed algorithm can be seen in the simulation results.
This paper considers the STAR-RIS in an ideal state, where the transmission and reflection phase-shift coefficients can be adjusted independently. However, in real applications, coupled transmission and reflection phase-shift coefficients are more common, and the related research problems are more complex. Therefore, the design of coupled STAR-RIS phase-shift coefficients is a crucial research direction for the future.